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Abstract 

In quantum cosmology the closed universe can spontaneously nucleate out of the state with 
no classical space and time. The semiclassical tunneling nucleation probability can be estimated 
as P ^ exp(— a 2 /A) where a=const and A is the cosmological constant. 

In classical cosmology with varying speed of light c(t) it is possible to solve the horizon 
problem, the flatness problem and the A-problem if c = sa n with s=const and n < —2. We show 
that in VSL quantum cosmology with n < — 2 the semiclassical tunneling nucleation probability 
is P ~ exp(— f3 2 A k ) with /3=const and k > 0. Thus, the semiclassical tunneling nucleation 
probability in VSL quantum cosmology is very different from that in quantum cosmology with 
c=const. In particular, it can be strongly suppressed for large values of A. In addition, we 
propose two instantons that describe the nucleation of closed universes in VSL models. These 
solutions are akin to the Hawking- Turok instanton in sense of 0(4) invariance but, unlike to 
it, are both non-singular. Moreover, using those solutions we can obtain the probability of 
nucleation which is suppressed for large value of A too. 

1 Introduction 

One of the major requests concerning the quantum cosmology is a reasonable specification of initial 
conditions in early universe, that is in close vicinity of the Big Bang. The three wave functions, 
describing the quantum cosmology has been proposed so far: the Hartle-Hawking's pQ, the Linde's 
[2], and the so-called tunneling wave function [3j. In the last case the universe can tunnel through 
the potential barrier to the regime of unbounded expansion. Following Vilenkin |4j lets consider the 
closed (k = +1) universe filled with radiation (w = 1/3) and A-term (w = —1). One of the Einstein's 
equations can be written as a law of a conservation of the (mechanical) energy: P 2 + U (a) = E, 
where P = —ad, a(t) is the scale factor, the "energy" E = const and the potential 

U(a) = c 2 a 2 (l-^ 
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where c is the speed of light. The maximum of the potential U(a) is located at a e = W3/2A where 
U{a e ) = 3c 2 /(4A). The tunneling probability in WKB approximation can be estimated as 

where d i < a>i are two turning points. The universe can start from a = singularity, expand to 
a maximum radius a£ and then tunnel through the potential barrier to the regime of unbounded 
expansion with the semiclassical tunneling probability (pQ). Choosing E = one gets a! i = and 
Oj = y3/A. The integral in ([I]) can be calculated. The result can be written as 

For the probability to be of reasonable value, for example P = 1/e ~ 0.368, one has to put A ~ 
0.3 x 10 65 cm -2 (see ([2])). In other words, the A-term must be large. On the other side, the universe 
once nucleated immediately begins a de Sitter inflationary expansion. Therefore the tunneling wave 
function results in inflation. And the A-term problem, which arises in this approach is usually being 
gotten rid of via the anthropic principle. In this case we have two Lorentzian regions (0 < a < a'^ 
a > ai) and one Euclidean region [a[ < a < a,{). The second turning point a = ai corresponds 
to the beginning of our universe. If A = then U(a) has the form of parabola and we get only 
one Lorentzian region. In this case, the universe can start at a = 0, expand to a maximum radius 
and recollapse. If E — > 0, the single Lorentzian region contracts to a single point, which lies in 
agreement with the tunnelling nucleation probability: P — > as A —* 0. However, as we'll show, in 
quantum cosmological VSL models the situation can be opposite, viz: the probability to find the 
finite universe short after it's tunneling through the potential barrier is P~ exp(— /3(n)A a ( n )) with 
a(n) > and (3{n) > when n < —2 or for —1 < n < —2/3. After the tunneling one gets the finite 
universe with "initial" value of scale factor ~ A -1 / 2 , so the probability to find the universe with 
large value of A and small value of is strongly suppressed. The reason for this lies in the behavior 
of potential U(a), which, for the case A — > 0, transforms into the hyperbola, located under the 
abscissa axis. As a result, such a universe can at a ~ start the regime of unbounded expansion. 
Therefore, we get the single Lorentzian region that doesn't contract to a point at E — > 0. 

This new property of VSL quantum cosmology will be discussed in next Section but new question 
arouse: the geometric interpretation of the quantum creation of a Universe with varying speed of 
light. We know that universe can be spontaneously created from nothing (when c = const) and this 
process can be described with the aid of the instantons solutions possessing 0(5) (if V((j>) has a 
stationary point at some nonzero value (f) = (fro = const) or 0(4) (as Hawking- Turok instanton [5]) 
invariance. So, what can be said about instantons in the VSL models? 

The whole plan of the paper looks as follows: in the next Section we'll consider the simplest 
VSL model: model of Albrecht-Magueijo-Barrow. Then we show that in framework of tunneling 
approach to quantum cosmology with VSL the semiclassical tunneling nucleation probability can 
be estimated as P ~ exp(— (3 2 A k ) with /3=const and k > 0. All corresponding calculations will 
be done for the case of the universe filled with radiation (w = 1/3) and vacuum energy. In the 
Section 3 we'll propose the non-singular instanton solutions possessing only 0(4) invariance (so 
the Euclidean region is a deformed four sphere). These solutions can in fact lead to inflation after 
the analytic continuation into the Lorentzian region. We will discuss these results in Sec. 4. 
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2 Albrecht-Magueijo-Barrow VSL model 



Lets start with the Friedmann and Raychaudhuri system of equations with k = +1 (we assume the 
G=const): 



d_ ^Gf^ + 3p^ + A(? _ 8nGp fc ^c\2 , Ac 2 



a o \ cr / o \ a / o v a 



a \" n 2 

c = eg I — = sa , p = wc p. 
\ a oJ 



3 ' 

(3) 



where a = a(t) is the expansion scale factor of the Friedmann metric, p is the fluid pressure, p is 
the fluid density, k is the curvature parameter (we put k = +1), A is the cosmological constant, cq 
is some fixed value of speed of light which corresponds to a fixed value of scale factor ao- 
Using (|3j) one gets 

3d / p\ cc(3 — a 2 A) 

Choosing w = 1/3 one can solve (JH) to receive 

_ M 3s 2 na 2 (™- 1 ) s 2 raAa 2n 
P ~ ^ + 87rG(n + !) " 87rG(n + 2 ) ' (5) 

where M > is a constant characterizing the amount of radiation. It is clear from the ([5]) that 
the flatness problem can be solved in a radiation-dominated early universe by an interval of VSL 
evolution if n < — 1, whereas the problem of A-term can be solved only if n < —2. The evolution 
equation for the scale factor a (the second equation in system ([3])) can be written as 

p 2 + U(a) = E, (6) 

where p = —aa is the momentum conjugate to a, E = 8irGA'I/3 and 

TT/ . s 2 a 2n + 2 2s 2 ka 2n+A 

U ^ = ^TT " W (7) 

The potential (JT]) has one maximum at a = a e = v / 3/(2A) such that 

s 2 3 n+1 

I7 e = C/(a e ) = 2 n+i A n+i (n + 1 ) (n + 2 )' (8) 

so C/ e > if (i) n < —2 or (ii) n > —1. The first case allows us to solve the flatness and "Lambda" 
problems. Another benefit of the model is a finite time region with accelerated expansion. 

2.1 The semiclassical tunneling probability in VSL models with n < —2: the 

case E <^U e 

One can choose n = —2 — m with m > 0. Such a substitution gives us the potential ([7|) in the form 

rr , s s 2 /2Aa 2 1 \ 
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Since © is similar to equation of movement of the particle of energy E in the potential ([9]), 
the universe in quantum cosmology can start at a ~ 0, expand to the maximum radius a- and then 
tunnel through the potential barrier to the regime of unbounded expansion with "initial" value 
a = di. The semiclassical tunneling probability can be estimated as 

P ~ exp ^-2 J 1 | p(a) \daj , (10) 

with 

I P(a) | = | p{a) |, | p{a) \= ^U m {a) - E, 

where E < U e . It is convenient to write E = U e sin 2 6, with < 6 < tt/2. 
For the case E <C £4 one can choose 



/ 3m fTf V^TT \ 1/m 

a l ~a 1 = J -r- — TT , a* ~ \ — — — — , (11) 



and evaluate the integral (|10|) as 



2(m + 1)A' V 2A V 

sin 



exp • (12 > 



where 



with 



I m (d) = I dzz~ 5 - 3m J- — , (13) 



3m m + 1 



z ^ d > = \ Ti — nT' ^ W = v L5 — — a — 

y 2(m + 1) \ sinfl 



l/m 



One can show that I m {6) > at < 6 <C 1. Thus, it is easy to see from ()12|) that the semiclassical 
tunneling probability P — > for large values of A > and P — > 1 at A — > 0. 

Note, that the case c=const can be obtained by substitution m = — 2 into the (I12p . Not 
surprisingly, this case will get us the well known result _P~ exp(— 1/A) (see [1]). 

2.2 The semiclassical tunneling probability with n < —2 and n > — 1 

In the case of general position the semiclassical tunneling probability with n = — 2 — m has the form 



/ s 3 A (3 m+ 4)/ 2 f^ dz^Ej^e) \ 

^ m ~exp^ 4 ^ 3(m+i)/2 ^___^ z3m+5 J, (14) 

where 

F m (z, 6) = -2 m+1 sin 2 6 z 2 ( m+1 ) + 2 x 3 m (m + l)z 2 - m3 m+1 , (15) 

z is dimensionless quantity and z[, z% are the turning points, i.e. two real positive solutions of the 
equation F m (z, 6) = for the given 6 (F m (z, 6) = does have two such solutions at < 8 < vr/2). 
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If m is the natural number then the expression (|14p has a more simple form. For example 

s 3 A 7 / 2 sin6> f Zi dz 
6irGhV2 J z > z i 



with 

V3 , V3 



2cos(0/2)' 4 2sin(0/2)" 
Similarly, P~ exp(— S 1 ), with 

s 3 A 5 sin6> f Zi dz 



where 



and so on. 

Therefore the probability to obtain (via quantum tunneling through the potential barrier) the 
universe in the regime of unbounded expansion is strongly suppressed for large values of A and small 
values of the initial scale factor aj = \/3/(2sin(0/2)\/A). In other words, overwhelming majority of 
universes born via the quantum tunneling through the potential barrier (|7|) have large initial scale 
factors and small values of A. 

Now, let us consider the case (ii), when n > — 1. The "quantum potential" has the form 

U(a) = s 2 a 2m (-- 2Aa2 - ) , (16) 
K ' \m 3(m + I) J ' v ; 

where m = n + 1 > 0. The points of intersection with the abscissa axis a are a>o = and a\ = 
\/3(m + l)/2Am. Choosing E = in equation ^ and substituting (|16|) into the (|10|) we get 



AirGh Jo W m 3(m + 1) 



with z\ = ^2>{m + l)/2m (The starting value z = means that the Universe tunneled from 
"nothing" to a closed universe of a finite radius ai = Zx/yA.). Thus, we have the same effect as if 
< m < 1/3. 

2.3 Peculiar cases with n = — 1 and n = — 2 

At last, lets consider the cases of n = — 1 and n = — 2. The formula (|14p is not valid in these cases 
(m = — 1 and m = 0) so we shall consider these models separately. 
If n = —1 (m = — 1) then 

_ M ; As 2 3s 2 a 
P ~^ + ^Ga* ~ ^Ga 1 g ^? 

therefore 

„ w ^( a ^(£)-^ +1 ), (17 , 
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where a* is constant and [a*]=cm. The potential (|17p has one maximum at a = a e = a/3/(2A) such 
that U e = U{a e ) = 2s 2 log(a e /a*), so if a e > a* then J7 e > 0. We choose a* = A _1//2 . This gives us 
U e = 0.41s 2 > 0. For the case E < (7 e the semiclassical tunneling nucleation probability is 



~ 6XP ("4^ /, W l0g * 2 - X + V ~ 6XP (-10^ J < (18) 

where the turning points are z\ = 0.721, Z{ = 1.812. As we can see from the f 1 1 8 f) . when n = —1 we 
receive the aforementioned effect again. 
If n = —2 (m = 0) then 

M s 2 K f a\ 3s 2 
P= -T + TT-Twlog — + 



a 4 2vrGo 4 V a */ 4vrGa 6 ' 
We choose a* = l/(a\/A), where a is a dimensionless quantity. Thus 

1 4A, / rr\ A 



17(a) = -s 2 ( ^ + ^ log (aaVAj + - ) . ( 19 ) 

The maximum of potential (|19p is located at the same point a e and 

u e=-—{ 3+los {— 

Therefore, U e > if a < 2e _3 / 4 /v / 6 ~ 0.386. Choosing a = 0.286 and E <C J7 e gets us the turning 
points 2$ ~ 0.77 and ~ 2.391. 

At last, the semiclassical tunneling nucleation probability is 



/ s 3 A 2 f Zi dz /l T " " f\ / 0.084s 3 A 2 

P ~ exp -— — — / r \l — k - - log{az) - - ~ exp 



AttGH L, 2 4 V z 2 3 bV ' 3 ^ 7rG/i 



3 Inst ant ons 



If we are going to describe the quantum nucleation of universe we should find the instanton solutions, 
simply putted as a stationary points of the Euclidean action. The instantons give a dominant 
contribution to the Euclidean path integral, and that is the reason of our interest in them. 
First at all, lets consider the 0(4)-invariant Euclidean spacetime with the metric 

ds 2 = c 2 (r)dT 2 + a 2 (r) (dip 2 + sin 2 ^ft 2 ,) . (20) 

In the case c = const one can construct the simple instantons, which are the 0(5) invariant four- 
spheres. Then one can introduce the scalar field </>, whose (constant) value 4> = (j)Q is chosen as 
the one providing the extremum of potential V(4>). The scale factor will be a(r) = H^ 1 sin Hr 
and after the analytic continuation into the Lorentzian region one will get the de Sitter space or 
inflation. Many other examples of non-singular and singular instantons were presented in [6] 

Now, lets consider the VSL model with scalar field. The corresponding Euclidean equations are: 



„ a' ,_ c?V' c 5 c'(Aa 2 -3) 2<f>'d 2cVc' 
+ a ~ ~1T + 4vrGaV + ~c <p~" 

a'\ 2 8itG fd) 12 2 \ c 2 Ac 2 
c 2 V + 



(21) 



3c 4 V 2 a 2 3 ' 
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where primes denote derivatives with respect to r. 

At the next step we represent the potential V in factorized form 

V = F(a)U((f>). (22) 

Indeed, lets for example consider the power- low potential ~ 4> k . If the coupling A is dimensionless 
one then we get 

n 



Since c = sa n then in the simplest case we come to 

Let 4> = 00 = const be solution of the (f2Tj) . (Note, that we don't require 4>q to be an extremum 
of potential.) Using the first equation of system (|2Tj) and (|22|) we get the equation for the F(a), 

_ 2 ^ F{a) = J^l a 4n-3 _ J^^-l (23) 

da a V ' AnGU 4ttGU ' V ' 

where Uq = U(4>o) = const. The integration of the ([23]) results in 



where C is the constant of integration and by assumption n/-l and n^O. Substitution of ([24 
into the second equation of the system (I2ip transforms it into the the model of nonlinear oscillator, 
integration of which result in 

a' 2 

— + u(a) = 0, (25) 

where 

. , uj 2 a 2 s 2 a 2n , , 

"<"' = — " W^Y <26> 

with uj 2 = 8irGUoC/(3s 2 ) and with the choice C > made. We can see that for c = const (i.e. 
n = 0) (|26p turns out to be an equation of the harmonic oscillator and we come to the well-known 
0(5) solution (but in this case </>o must be the stationary point of V). 

Equation (|25p naturally describes the "movement of a classical particle" with zero-point energy 
in mechanical potential (|26|) , Depending on value of n this potential can take one of four distinct 
forms (excluding the well-known classical case n = 0, which lies beyond the scoop of this article). 

Case 1: n < 0. Here we have one Euclidean (0 < a < a\) and one Lorentzian (a > oi) regions 
where 

«i=(— «=) • ( 27 ) 



On the bound between Euclidean and Lorentzian regions {a = a\) we have a 1 = 0. 

This mechanical potential is unbounded from below at a — > 0. With this in mind, we'll have to 
ascertain that the Euclidean action for our solution will stay finite. The gravitation action has the 
form 

Sgrav = - J d A x£^y/gR. 
We are using the dimensionless variables x° = cor/ao, x 1 = ip and so on. Calculating R we get 



R 



c 2 a 2 



c - (^-J ((l-n)a +aa) 



(28) 
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so we do have the potential divergence at a = 0. Multiplying (|28p on the yfg and c 3 and using the 
equation of motion we get the expression: 



/ a 2n+3 
R^c 3 ~ 6c (2 - n)Lo 2 ^—^ - - 

V a 1 



2n+3 „„2 „4n+l 



4n — 1 / ' 

nap ' 



(29) 



where the most dangerous multiplier factor is a 1+4n . But if— 1/4 < n < then the Euclidean 
action becomes finite and therefore, we end up with the legitimate gravitation instanton. In a 
similar manner, using (|22p and f|24j) we get for the scalar field (in dimensionless x^): 



.1-371 



C Q /o, ,2„3(l+n) , 6rlC a _ fi£0 



3nc2a 1+5n Ac2a 3+5n , 



8vrG v (n-l)ao 



2n „2n 







therefore the instanton exists for n > —1/5. This requirement is stronger than the one for the 
gravitation instanton where n > —1/4 (see (|29p). 

Case 2. < n < 1. Here the potential u(a) suffers no singularity at a = 0, but u(0) = 0. Also 
this potential has a minimum at 



a 



, 1/(1-™) 

n 



lo V 1 — n 



and is equal to zero at (I27h . hence, once again creating one Euclidean and one Lorentzian regions, 
separated by (j2T|) . 

Case 3. n = 1. This case is somehow special, since for such n the solution of (j23[) shall be 



instead of ([24]) . and hence, the equation of (|26j) shall be substituted by 

u(a) = a 2 - s 2 In aj . (30) 

2 

It is easy to see that this function has two zeros (at a\ = and ci2 = exp(^y)), is strictly positive on 
interval (a%, 02) and strictly negative outside of it. Therefore, this case doesn't allow an instanton. 

Case 4. n > 1. The potential u(a) is strictly positive. The instanton doesn't exist either. 

Both of a newly founded solutions possess only 0(4) invariance just like Hawking- Turok instan- 
ton (so the Euclidean region is a deformed four sphere) but, unlike to it, they are all non-singular. 
Note that if the value a is sufficiently large then one can neglect the second term in (I26|) (after 
the analytic continuation into the Lorentzian region) therefore, as in the case of the usual 0(5) 
instanton, one can get the de Sitter universe, i.e. the inflation. 

The equation ([25|) has no terms with A. In other words, the scale factor a(r) doesn't depend 
on the value A (although being dependant on the Uq). Therefore, the full Euclidean action S E = 
Sgrav + Afield nas the- form, 

S E = So — ASi, 

where Sq and S\ are both independent of the A. Returning to what has been said in Introduction, 
there exist three common ways to describe the quantum cosmology: the Hartle-Hawking wave 
function exp(—S E /h), the Linde wave function exp(+5' s /^) and the tunneling wave function. In 
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the second Section we have been working with the tunneling wave function. In case of instantons 
situation becomes slightly different. If Si > then (as a first, tree semiclassical approximation) 
we should choose the Linde wave function, whereas for the case S\ < the Hartle-Hawking wave 
function seems more naturally. 

In conclusion, we note that another choice of C (C < and C = 0) eliminates all possible 
instantons. 

4 Discussion 

VSL models contain both some of the promising positive features [7] and some shortcomings and 
unusual (unphysical?) features as well [8]. But, as we have shown, application of the VSL principle 
to the quantum cosmology indeed results in amazing previously unexpected observations. The first 
observation is that the semiclassical tunneling nucleation probability in VSL quantum cosmology is 
quit different from the one in quantum cosmology with c=const. In the first case this probability 
can be strongly suppressed for large values of A whereas in the second case it is strongly suppressed 
for small values of A. This is interesting, although we still can't say that VSL quantum cosmology 
definitely results in solution of the A-mystery. The problem here is the validity the WKB wave 
function. And what is more, throughout the calculations we have been omitting all preexponential 
factors (or one loop quantum correction) which can be essential ones near the turning points. 
Another troublesome question is the effective potentials in VSL models, being unbounded from 
below at a — > 0. The naive way to solve this problem is to use the Heisenberg uncertainty relation 
to find those potentials with the ground state. However, this is just a crude estimation. To describe 
the quantum nucleation of universe we have to find the instanton solution which, being a stationary 
point of the Euclidean action, gives the dominant contribution to the Euclidean path integral. As 
we have seen, such solutions indeed exist in VSL models. Those instantons are 0(4) invariant, are 
non-singular, and provide an inflation as well. They describe the quantum nucleation of universe 
from "nothing" and, what is more, upon usage of these solutions we can obtain the probability of 
a nucleation which is suppressed for large value of A using either Linde or Hartle-Hawking wave 
function. 

Note, that we can weaken the condition n > —1/5 to obtain a singular instanton suffering 
the integrable singularity (i.e. such that the instanton action will be finite) in the way of the 
Hawking- Turok instanton. However, there exist some arguments [9], that such singularities, even 
being integrable, still lead to serious problems with solutions. 

In conclusion, we note that obtained instantons both have a free parameter (lo 2 ) so we are free 
to use the anthropic approach to find the most probable values of A too. 

Acknowledgements: After finishing this work, we learned that T.Harko, H.Q.Lu, M.K.Mak and 
K.S.Cheng [10], have independently considered the VSL tunneling probability in both Vilenkin and 
Hartle-Hawking approaches. The interesting conclusion of their work is that at zero scale factor the 
classical singularity is no longer isolated from the Universe by the quantum potential but instead 
classical evolution can start from arbitrarily small size. In contrast to |10j . we attract attention to 
the problem of A-term and instantons in VSL quantum cosmology. 
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